Ultrafast spatio-temporal dynamics of terahertz generation by ionizing two-color 

femtosecond pulses in gases 
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We present a combined theoretical and experimental study of spatio-temporal propagation effects 
in terahertz (THz) generation in gases using two-color ionizing laser pulses. The observed strong 
broadening of the THz spectra with increasing gas pressure reveals the prominent role of spatio- 
temporal reshaping and of a plasma-induced blue-shift of the pump pulses in the generation process. 
Results obtained from (3+l)-dimensional simulations are in good agreement with experimental 
findings and clarify the mechanisms responsible for THz emission. 
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Far-infrared radiation in the THz range has devel- 
oped into a sensitive probe of low-frequency excitations 
of condensed matter and into an analytical and imag- 
ing tool for a broad range of applications. Recently, in- 
tense THz pulses have been applied for inducing nonlin- 
ear light-matter interactions and for studying quantum- 
coherent charge transport phenomena in solids. In this 
new area of research, the generation of well-defined THz 
field transients with a high amplitude represents a key 
issue. Both electron accelerator- and laser-based sources 
have been developed to generate THz transients. Con- 
ventional laser-driven THz sources are based on semicon- 
ductor photoconductive switches [H and nonlinear fre- 
quency conversion in crystals , providing a comparably 
small THz field strength in a spectral range limited by the 
absorption of the crystals. In an alternative approach Q, 
a short pump pulse at 800 nm and its second harmonic at 
400 nm are focused into a gaseous medium to generate a 
plasma. Intense THz pulses with field amplitudes as high 
as 400 kV/cm and remarkably broad spectra have been 
reported with this method [3H13|. However, the physi- 
cal mechanisms behind the observed THz generation re- 
main controversial. Initially, THz generation has been 
explained by four- wave mixing rectification 0, [l4| and 
experimental results have been interpreted along those 
lines (see, e.g., 0-0, H, EH)- In contrast, other studies 
attribute THz emission to a laser- induced plasma current 
in the asymmetric two-color laser field 0, E3, So far, 
theoretical work has mainly focused on the analysis of lo- 
cal fields by using photocurrent models. An alternative 
approach is based on particle in cell simulations [l5r{13] 
which are computationally very expensive and limited to 
propagation lengths of a few micrometers only. 

In this Letter we present a combined theoretical and 
experimental study of THz generation by ionizing two- 
color femtosecond pulses in a gas. Extensive numeri- 
cal simulations were performed using for the first time a 



(3+l)-dimensional code based on an unidirectional pulse 
propagation equation which includes the plasma dynam- 
ics responsible for the observed THz generation. This 
approach allows a comprehensive description of the prop- 
agation of all fields taking into account their spatio- 
temporal reshaping induced by the plasma effects and the 
optical Kerr nonlinearity. We demonstrate that spatio- 
temporal propagation effects are indispensable for under- 
standing the generation process and influence the THz 
spectrum substantially. In both experiments and simu- 
lations we observe a remarkable broadening of the THz 
spectrum with increasing gas pressure. Such broadening 
is a result of a sensitive dependence of the THz spectrum 
on small phase and frequency shifts induced by nonlin- 
ear propagation of the fundamental and second harmonic 
pulses. We show that the stepwise character of the ion- 
ization process on the sub-femtosecond time scale is es- 
sential for such dependence. 

Our THz plasma source (inset of Fig. [1]) is driven by 
40 fs pulses (800 nm) with pulse energies of ~ 300 //J at 
a repetition rate of I kHz from a Thsapphire laser sys- 
tem. Intensities far above the field ionization threshold 
are reached by focusing the beam of 8 mm diameter with 
an achromatic lens (L) of 40 mm focal length. A 0.1 mm 
thin /3-barium borate (BBO) crystal (C) cut for type I 
second-harmonic generation is additionally inserted into 
the convergent beam 7 mm before the focus. The setup 
is placed in a closed chamber filled with argon at various 
pressures between 1 and 1000 mbar. THz radiation emit- 
ted from the plasma volume in the focal spot is collected 
by a parabolic mirror (M) with a diameter of 25.4 mm 
at a distance of its effective focal length of 12.7 mm. 
The THz spectrum is measured by a Michelson interfer- 
ometer. Intensity interferograms were recorded with a 
mercury cadmium telluride (HgCdTe) detector by vary- 
ing the path difference between the two arms, and the 
THz spectra were obtained by Fourier transformation. 
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FIG. 1: (color online) The mechanism of THz generation: (a) 
The two-color electric field E generates free electrons with a 
step-wise modulation of the electron density p via tunneling 
photoionization. The ionization occurs mostly near the max- 
ima of electric field at time instances t n . (b) This leads to 
a slow component of the current J growing on time scales of 
the pulse duration, which acts as a source for THz emission. 
Inset in (a): scheme of the experimental setup. 



For the description of the nonlinear propagation of 
electromagnetic fields valid in the full range from THz 
down to UV/VUV wavelengths beyond paraxial and 
slowly-varying envelope approximations we use the fol- 
lowing unidirectional pulse propagation equation [2l| for 
linearly polarized pulses 
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is the Fourier transform (indicated 



Here, E(k x ,k y ,z 
by " ) of the electric field with respect to x, y, and t, 
k(u) = ujn(u))/c is the wavenumber, lo the frequency, c 
is the speed of light and n(oj) is the refractive index of 
the gaseous medium, in our case argon 19]. The first 
term on the r.h.s. of Eq. (fT]) describes linear dispersion 
and diffraction of the pulse. The nonlinear polarization 
"Pnl = -Fkcrr + iJ e /uj + iJi oss /uj originates from the op- 
tical Kerr effect Pkctt, the electron current J e and a loss 
term Ji oss due to photon absorption during ionization. 
In -Pkcit we take into account different nonlinear suscep- 
tibilities for neutral atoms and ions [2(|. However, the 
Kerr response has negligible influence on the THz spec- 
trum. The plasma dynamics is described by the electron 
density p e (t), obeying 



p e (t)=W ST (E)[p at -p e (t)}, 



(2) 



where p at denotes the neutral atomic density and dot 
the time-derivative. We use the quasi-static tunneling 



ionization rate for hydrogen-like atoms given by [10| 
W ST (E) = 4w a (r H ) 2 - 5 [£ a /|£|]exp[--2(rtf) 1 - 5 £j3|£|], 
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Uh and U Ar are the ionization potentials of hydrogen and 
argon; m e and q are the electron mass and charge, respec- 
tively. For the anticipated electric field strengths, we ne- 
glect multi-photon and avalanche ionization. The trans- 
verse macroscopic plasma current J e (t) is determined by 
the microscopic velocity distribution v(t,to) of electrons 
born at the time to 
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v(t,t )p e (ta)dt . 



(3) 



Assuming zero velocity for new-born electrons and ne- 
glecting the influence of the magnetic field and electron- 
electron interaction, the electron velocity reads v(t, to) — 
J t E(t) exp[— v e (t— r)]dr, where v e is the electron-ion 
collision rate. Substituting this into Eq. © we obtain 



J e {t) + v e J e {t) = —E{t)p e (t). 



(4) 



The above analysis demonstrates that we regain a well- 
known equation for the current (see, e.g., |20(). To for- 
mally ensure energy conservation during ionization, the 
additional term Ji oss = VFst [i?] (Pat — p e )V Ar/E is intro- 
duced in Eq. (Jl]) such that the energy dissipation J\ oss E 
equals the ionization energy loss. 

Let us now illustrate the mechanism behind THz gen- 
eration. In Fig. [IJa) we present the electron density (in 
green), which shows a stepwise increase near the tunnel 
ionization events at the field maxima (red curve). Such 
behavior has been recently observed in real-time exper- 
iments with sub- femtosecond resolution jl8| . In a sim- 
plified model, we can assume rectangular steps in the 
electron density [p e (t) = T, n p n 5(t — t n )] and v e = 0, and 
thereby obtain a discrete version of Eq. ([3j> : 



J e (t) ~^p n H{t - t n )[v f (t) ~ v f (t n )], 



(5) 



where H(t) is the Heaviside step function, Vf(t) = 
— f E(T)d,T is the free electron velocity such that 
v(t,t n ) = [vf(t) — Vf(t n )}, and p n is the electron density 
created in the nth ionization event. For a monochromatic 
electric field Vf(t n ) = 0. The Fourier transformation of 
the step function H(t) exhibits a low- frequency spectrum 
~ therefore THz radiation is generated by the terms 
in Eq. §5§ proportional to Vf(t n ), while the terms pro- 
portional to Vf(t) contribute in the spectral range of the 
pump fields. 

In order to model the experimental conditions in the 
plasma spot, we consider a Gaussian input beam with a 
waist wo = 100 /xm and a duration of the 800 nm pump 
pulse of (FWHM) t p = 40 fs (pulse energy 300 pJ). The 
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FIG. 2: (color online) (a) Iso-electron-density surface at p e — 
5 x 10 17 cm -3 for 200 mbar gas pressure. Computed THz 
fields (u> < 500 ps -1 ) at (b) z — 0.2 mm and (c) z = 1 mm. 



energy of the second harmonic at 400 nm is chosen as 
12% of the fundamental as estimated from the experi- 
ment. The duration t p and the waist wq of the second 
harmonic are by a factor \/2 smaller than the values for 
the fundamental. The phase angle between the two com- 
ponents is zero initially and shifts during propagation 
to non-zero values. The pulsed input beams are focused 
(/ = 0.5 mm, in order to have a comparable ratio wq/ f in 
experiment and simulation) into the argon atmosphere. 
Figure [2] shows (a) the evolution of the computed plasma 
density and the transverse distribution of the THz field at 
(b) 0.2 mm and (c) 1 mm for a gas pressure of 200 mbar. 
The field intensities reach the ionization threshold shortly 
after the starting point of the simulation at z — mm, 
and form a 0.7 mm long plasma channel. THz fields in- 
side this focal region reach values of the order of GV/m 
and exhibit a strong diffraction. 

Let us now investigate the dependency of THz gen- 
eration on the gas pressure. Measured spectra over the 
complete pressure range between 1 mbar and 1000 mbar 
are shown as a contour plot in Fig. EJa). Our HgCdTe 
detector is sensitive in a frequency range uj = 2ttv from 
125 ps -1 to 1070 ps -1 , i.e., the drop of the measured 
energy spectral density around 125 ps -1 is due to the 
decreasing sensitivity. For a comparison of theory and ex- 
periment, the high-frequency part of the spectra is most 
relevant. In our setup, we mainly detect THz radiation 
generated in the focus of the mirror (M). Simple ray trac- 
ing estimates indicate that the length of this focal region 
is < 0.3 mm. The almost vanishing spectrum at small 
argon pressure shows clearly that the plasma and not the 
BBO crystal acts as a source of the emitted radiation. In 
the region from zero to ~ 300 mbar the spectral width 
increases strongly. The highest frequencies even beyond 
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FIG. 3: (color online). Measured THz spectra (a) and simula- 
tion results (b) for pressures between 1 and 1000 mbar. In (c), 
experimental (solid lines) and theoretical (dashed lines) spec- 
tra are compared for 100, 200 and 400 mbar (green, blue and 
red curves, resp.). In (b),(c), the part below uj — 125 ps -1 is 
shaded, because no experimental values are available. In (d), 
the overall THz yield versus pressure is shown (dashed green 
line represents the simulation, solid line the experiment). 



300 ps -1 are detected at pressures larger than 300 mbar. 
Above ~ 500 mbar the slope of the high frequency wing 
stays rather constant. The measured THz yield [solid 
line in Fig. [3][d)] grows linearly with the pressure up to 
400 mbar before it saturates. 

In Fig. [3][b,c), THz spectra computed with our simula- 
tion code are shown. The simulated spectra in Fig.[3][b,c) 
are obtained by integration over the transverse coordi- 
nates (x,y). We find very good agreement between ex- 
periment and simulation below 500 mbar for THz fields 
generated at the beginning of the plasma spot around 
z = 0.2 mm (0.3 mm before the linear focus). THz fields 
generated upon further propagation become spectrally 
much broader. Hence, we conjecture that the parabolic 
mirror in the experiment images the leading part of the 
plasma spot only. The calculated THz yield increases 
linearly with gas pressure [dashed line in Fig. [Ud)], in 
agreement with experimental results up to 400 mbar. 
The saturation of the experimental yield at higher pres- 
sure is likely due to additional THz losses upon further 
propagation towards the mirror, whereas the simulated 
yield is computed directly at the position z = 0.2 mm. 

The observed pressure dependence of the THz spectral 
maximum and width gives insight into important features 
of plasma-induced THz generation. The dependence of 
the spectral width on pressure can not be explained by 
the local plasma current in which the variation of pres- 
sure results only into an amplitude scaling of the current. 
Instead, it originates from pressure dependent nonlinear 
propagation effects. For the intensity range and plasma 
interaction length of the experiment, the calculated spec- 
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FIG. 4: (color online) (a) On-axis (x = 0,y — 0) spectral 
intensity I(to) and (b,c) temporal shape of E(t) for 400 mbar 
gas pressure. In (a), I(oj) at z = 0.1 mm and z = 0.2 mm 
are shown by the solid red and dashed blue line, respectively. 
The temporal profiles of the field E(t) (blue curve) and its 
low frequency part (red curve) are shown in (b) for z — 0.1 
mm and in (c) for z = 0.2 mm. I{uj) and E(t) at z = 0.1 mm 
[red curves in (a) and (b)] are amplified by a factor 100 and 
10, respectively, in order to improve visibility. 



propagation distance can be explained by propagation 
effects of the pump fields. We would like to stress that 
phase and frequency relation of fundamental and second 
harmonic fields always change during propagation. 

In conclusion, we presented a theoretical and experi- 
mental investigation of photo-current induced THz gen- 
eration. Using (3+l)-dimensional simulations and exper- 
imental measurements we show a strong dependence of 
the THz spectra on gas pressure. Our results give insight 
into the important influence of nonlinear propagation ef- 
fects and the mechanism of THz generation. Plasma- 
induced blue-shifts of the driving pulses play a key role 
in the pressure dependent broadening of the THz spec- 
tra and confirm that the emission process of THz radia- 
tion is associated with a stepwise modulation of the tun- 
neling ionization current. We believe that our findings 
open interesting perspectives to control THz emission in 
a broader spectral range. 

This work has been performed using HPC resources 
from GENCI-CINES (Grant 2009-x2009106003). 



tral evolution of the pump pulses at 400 nm (u = 2tt x 750 
THz) and 800 nm (lo = 2tt x 375 THz) shows that their 
spectral broadening is negligible. We observe, however, 
small blue-shifts 5lu of the central frequencies caused by 
the nonlinear plasma-induced change of the refraction in- 
dex [22|, [23|. These shifts are ~ 6 ps _1 in the fundamental 
and ~ 2.5 ps _1 in the second harmonic at z = 0.2 mm 
for 400 mbar, and depend strongly on the gas pressure. 

Surprisingly, these very small frequency shifts have a 
dramatic influence on the generated THz spectrum. This 
effect is most pronounced in the on-axis spectra, where 
the intensity is maximal. In Fig. QJa), the calculated 
on-axis THz spectra are plotted for 400 mbar gas pres- 
sure, displaying dramatic changes in spectral shape and 
the maximum frequency between 0.1 mm (red curve) and 
0.2 mm (blue curve). Figs. 4 (b,c) show the correspond- 
ing time-dependent field and its low- frequency part. To 
explain the physical origin of this very sensitive depen- 
dence, let us go back to our simplified model [Eq. ((SJ)]. 
For a two-color optical field, E = A\ cos[(wo + 8uj)t] + 
A2 cos(2woi + 0), the field maxima in every half-cycle are 
given by u>ot n nir — nir5u>/uJo — (— l) n 2r sin^, provided 
that A2/A1 = r <C 1 and n5cj <C uiq. Hence, the points 
in time t n and therefore the free velocities Vf(t n ) alter 
significantly when Sto (and 9) change upon the propa- 
gation. As seen above, the low-frequency spectrum is 
determined by a sum over contributions ~ Vf(t n ), and 
this sum finally determines the THz spectral shape in 
Fig. Eta). In the full (3+l)-dimensional geometry, the 
spectral shapes generated at different spatial points are 
added and averaged, leading to the strong spectral broad- 
ening observed in Fig. Etc). Thus, the above described 
dependence of the THz spectral shapes on pressure and 



[1 

[2 

[3; 

K 
[5 
[6 
[T 
[8 

[9 

[10 

[11 

[12 

[13 

[14 
[15 

[16 

[17 
[18 
[19 

[20 

[21 

[22 

[23 



D. You, R. R. Jones, P. H. Bucksbaum, and D. R. Dykaar, 
Opt. Lett. 18, 290 (1993). 

D. H. Auston, K. P. Cheung, J. A. Valdmanis, and D. A. 

Kleinman, Phys. Rev. Lett. 53, 1555 (1984). 

D. J. Cook and R. M. Hochstrasser Opt. Lett. 25, 1210 

(2000). 

P. Zhang et al., Phys. Rev. Lett. 91, 225001 (2003). 

M. Kress et al. Opt. Lett. 29, 1120 (2004). 

T. Bartel et al., Opt. Lett. 30, 2805 (2005). 

M. Kress et al., Nat. Phys. 2, 327 (2006). 

Xu Xie, J. Dai, and X.-C. Zhang, Phys. Rev. Lett. 96, 

075005 (2006). 

K. Kim, J. Glownia, A. Taylor, and G. Rodriguez, Opt. 
Expr. 15, 4577 (2007). 

M. D. Thomson, M. Kress, T. Loefner, H. G. Roskos, 

Laser & Photon. Rev. 1, 349 (2007). 

K. Reimann, Rep. Progr. Phys. 70, 1597 (2007). 

K. Kim, J. Glownia, A. Taylor, and G. Rodriguez, Nat. 

Phot., 2, 605 (2008). 

A. Houard, Yi Liu, B. Prade and A. Mysyrowicz, Opt. 
Lett. 33, 1195 (2008). 

J. Penano et al, Phys. Rev. E 81 026407 (2010). 
H.-Ch. Wu, J. Meyer-ter-Vehn, and Zh.-M. Sheng, New 
J. Phys. 10, 043001 (2008). 

M. Chen, A. Pukhov, X.-Yu Peng, and O. Will, Phys. 
Rev. E 78, 046406 (2008). 

W.-M. Wang et al., Opt. Expr. 16, 16999 (2008). 
M. Uiberacker et al., Nature 446, 627 (2007) . 
A. Dalgarno and A. E. Kingston, Proc. Royal Soc. Lon- 
don A 259, 424 (1960). 

P. Sprangle, J. R. Penano, B. Hafizi, and C. A. 

Kapetanakos, Phys. Rev. E 69, 066415 (2004). 

M. Kolesik and J. V. Moloney, Phys. Rev. E 70, 036604 

(2004). 

W. M. Wood, C. W. Siders, and M. C. Downer, Phys. 
Rev. Lett. 67, 3523 (1991). 

S.C. Rae and K. Burnett, Phys. Rev. A 46, 1084 (1992). 



